Zero- Temperature Theory of Collisionless Rapid Adiabatic Passage from a 
Fermi Degenerate Gas of Atoms to a Bose-Einstein Condensate of Molecules 
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We theoretically examine a zero-temperature system of Fermi degenerate atoms coupled to bosonic 
molecules via collisionless rapid adiabatic passage across a Feshbach resonance, focusing on satu- 
ration of the molecular conversion efficiency at the slowest magnetic-field sweep rates. Borrowing 
a novel many-fermion Fock-state theory, we find that a proper model of the magnetic-field sweep 
can systematically remove saturation. We also debunk the common misconception that many-body 
effects are responsible for molecules existing above the two-body threshold. 
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/n£roduciio?i.-Magnetoassociation creates a molecule 
from a pair of colliding atoms when one of the atoms 
spin flips in the presence of a magnetic field tuned near 
a Feshbach resonance . Recently, ultracold 0, and 
condensate molecules have been created via magne- 
toassociation of a Fermi gas of atoms, in the course of 
efforts to create superfiuid Cooper-paired atoms @, ||J 
(see also Refs. 0). The backbone of these experiments is 
rapid adiabatic passage: the ground state of the Feshbach 
system is all atoms far above the molecular-dissociation 
threshold and all molecules far below it, so that a slow 
sweep of the magnetic field from one extreme to the other 
converts atoms into diatomic molecules. 

Finite-temperature mean-field theory of magnetoasso- 
ciation of a Fermi gas of atoms leads to two types of insta- 
bilities against molecule formation. One is the thermo- 
dynamic instability of a Fermi sea against the formation 
of Cooper pairs § , a trait of superconductors whose ana- 
log is passed on to Feshbach-resonant superfluids |9(. A 
thermodynamical instability occurs because pairing low- 
ers the energy, and coupling to a reservoir with a low 
enough temperature leaves the system prone to pairing. 
The other is a dynamical instability, whereby the larger 
state space of the molecules, owing somewhat to Pauli 
blocking, leaves the atoms prone to spontaneous associa- 
tion [lOj . The role that temperature plays in this process 
is an open question experimentally, as well as a matter 
of theoretical contention. 

Physically ^(j > high temperature lessens the chance of 
an atom occupying an arbitrary level in the Fermi sea, the 
dynamical instability becomes less effective and the effi- 
ciency of even the slowest rapid adiabatic passage there- 
fore saturates (c.f., Fig. QJ. The mean-field theory be- 
hind this understanding agrees semi-quantitatively with 
experiments 0]; nevertheless, a recent zero-temperature 
Landau-Zener theory predicts that saturation is funda- 
mental to the collisionless regime If temperature is 
not a limiting factor, then any zero-temperature model 
of collisionless rapid adiabatic passage should ultimately 
display saturation, e.g., a Fock-state approach similar to 
the theory of cooperative association of Bose- Fermi mix- 
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FIG. I: Predicted temperature dependence for the efficiency 
(|/3| 2 ) of collisionless rapid adiabatic passage from quantum 
degenerate 40 K atoms to 40 K2 molecules. As the temperature 
in Fermi units (fieF/fes) is lowered, atoms become more likely 
to be affected by the dynamical stability that forms molecules. 
The magnetic field was swept linearly at the (inverse) rate 
1/B = 400^s/G. Figure reproduced from Ref. [fu- 



tures of atoms into Fermi molecules [l^ . Unfortunately, 
computing power is presently sufficient for calculations 
with only about 20 atoms total at best, precluding any 
brute-force 01 test of saturation. Here we apply a novel 
large- fermion-number theory |13| to demonstrate near- 
unit-efficient collisionless rapid adiabatic passage in the 
limit of zero temperature, thereby ruling out any funda- 
mental ceiling to the molecular conversion, and bolstering 
Ref. [13 (and also Ref. 0). 

This development is outlined as follows. After briefly 
introducing the collisionless model, we focus on rapid 
adiabatic passage and confirm the reduced-space map- 
ping |l3j by comparison with exact few-particle results. 
Increasing the total particle number to 2 x I0 2 , we then 
observe what appears to be saturation at about ~ 50%. 
However, including fluctuation effects in the rate at which 
the system is swept across the Feshbach resonance, we 
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find that saturation can be systematically removed, and 
near-unit efficiency can be achieved for any particle num- 
ber. Lastly, from the single pair results we also debunk 
the commonly held notion that many-body effects are re- 
sponsible for the existence of molecules above the thresh- 
old for molecular dissociation. 

Collisionless Gas Model.-We model an ideal two- 
component gas of fcrmionic atoms coupled by a Fesh- 
bach resonance to bosonic molecules. In the language 
of second-quantization, an atom of mass m and momen- 
tum Kk is described by the annihilation operator 0^,1(2), 
and a molecule of mass 2m and similar momentum is 
described by the annihilation operator b^. All operators 
obey their (anti)commutation relations. The microscopic 
Hamiltonian for such a freely-ideal system is written 



Hamiltonian J5J yields 

iC m = [N-m]6C r . 



VN-m + lD™- 1 ^-! 
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Here C m (t) = C 

N—m,rii , . ...tln 

(t) is a column vector of 
all the amplitudes corresponding to the (^) possible 
arrangements of m atom pairs among the N available 
fermion modes, and Dj is an ( I j)x( 1 j) dimensional ma- 
trix that contains only unit and zero elements determined 
by Ci and Cj. The problem with the system Q is that 
there are 2 N amplitudes, which limits most numerical 
experiments in rapid adiabatic passage to about N = 10 
(see also Refs. however, by multiplying Eqs. Q by 

the appropriate column vector u mj Ar, any redundant am- 
plitudes can be eliminated [T^- The remaining N + 1 
amplitudes evolve in time according to [T^ | 



where repeated greek indices imply a summation (a = 
1,2). The free-particle energy is Hek = H 2 k 2 /2m, the 
atom (molecule) chemical potential is fi./V( mo i), an d the 
detuning S is a measure of the binding energy of the 
molecule (8 > is taken as above threshold), the mode- 
independent atom-molecule coupling is n cx 1/VV with 
V is the quantization volume. 

We have already imposed the ideal conditions for atom- 
molecule conversion with /xi = /i2 = /i. An appropriate 
unitary transformation then shuffles [i into the defini- 
tion of /Xmoi which, in turn, can be absorbed into the 
detuning and written off as an effectively dc bias (see 
also Refs. ^3)- Since magnetoassociation usually occurs 
much faster than any trapped-particle motion, an explicit 
trap can be neglected along with the free-particle ener- 
gies €k lla. For the sake of simplicity, and to compare 
with Ref. [T(| , we neglect all molecular modes except the 
k + k' = mode, bo = b, so that 



= [N - m]6a m 

+ K [y/m(N- 



1 + l)a m _i 
+y/m + 1 (N - m)a m +i] 



(5) 
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= 5tfb + K^2 [& + ak,ia- k ,2 + H.c] 



Absent losses, the total particle number is conserved, 
2(&t&) + X)k( a L a k^> = 2n + ^ kCT n k<r = 2N, where n 
is the number of molecules, n^.a = 0, 1 is the number of 
atoms per mode (k) per species (a), and 2N is the total 
number of atoms were all the molecules to dissociate. 
For a fixed number of particles equal to the number of 
fermion modes, the Fock-state wavefunction is fl3| 

N 

IV>W) = C N~m>, ni ,...,n N 
m'=0 {n fe } 

(3) 

The time dependence of the system is determined by 
the Schrodinger equation, ihdt^) = H\ip), so that the 



where the sum of all (^ ) amplitudes with N — m 
molecules and m free atom pairs is defined as yj (^ ) a m = 
u m ,ArC m = J2n k CN-n hni ,...,n N (with a m normalized to 
the number of permutations of m atoms in N states). 
Lastly we will need the molecular fraction |/3| 2 = 

2<6t6)/(2JV) = (1/JV) EZ=o(N ~ m)\a m \ 2 . 

Rapid Adiabatic Passage-Putting fluctuations pL7| 
momentarily aside, the relevant frequency scale is Q = 
VNk oc y/p [H E3, where yfp is the so-called collec- 
tive enhancement factor. "Adiabatic " is therefore defined 
qualitatively as the detuning changing by an amount Jl 
in a time I/O, or |<j| < fi 2 . Modeling the time dependent 
detuning as 6 = — ££l 2 t, sweeps with £ ~ 1 should qualify 
as adiabatic. Off hand, Fig. Ola) confirms this intuition 
for N — 4; also, noting that the full results are shifted 
for clarity, the reduced system J5J indeed reproduces the 
full system (J2J. Making a more full use of the reduced- 
space theory, Fig. Gib) illustrates that the efficiency of 
rapid adiabatic passage in fact decreases for increasing 
particle number, saturating at about 50% for N = 10 2 . 
Nevertheless, if we account for fluctuations, then the rel- 
evant frequency scale is £1/ In N 17]. Now the detuning 
should change by SI/ In TV in a time (Sl/lniV)" 1 , sug- 
gesting the detuning-sweep model S(t) — — £(f2/ In N) 2 t. 
Indeed, Fig.^c) shows that the N = 10 2 and N = 1 re- 
sults agree nicely, and are absent any evident saturation. 

We can also make a rough comparison with the zero- 
temperature limit in Fig.^ Magnetic fields are converted 
(t).\N — m' ,ni, . . . ,n N ) into detunings according to S = A^(B — B )/h, where 

the difference in magnetic moments between the atom 
pair and a molecule is A M , and Bq is the magnetic-field 
position of resonance. For N = 10 5 atoms of 40 K in a 
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typical [2\ trap the peak density 
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FIG. 2: Molecular condensate fraction as a function of the 
detuning in rapid adiabatic passage across a Feshbach res- 
onance, beginning above threshold (6 > 0). (a) For N = 4, 
comparison of solution to the full equations of motion [Eq. ijljl. 
dashed line] with the solution to the reduced-space equations 
of motion [Eq. 10, solid line]. The full results are shifted 
for clarity; the two calculations are otherwise indistinguish- 
able. The detuning sweep model is S(t) = — £fi 2 t, with £ = 1. 
(b) Using the same sweep model, we find apparent saturation 
for increasing particle number: N = 4 (solid line), N = 10 
(dashed line), and A^ = 10 2 (dotted line), (c) A fluctuation- 
adjusted sweep model, S(t) — — £(fi/ In A") 2 t, leads to near- 
unit efficiency for N = 10 2 (solid line), as well as solid agree- 
ment with the N — 1 results. Here the dimensionless sweep 
rate is again unity, £ = 1. (d) Results for A^ = 10 2 and 
£ = 7.6, an estimate of a zero-temperature sweep for 40 K (see 
text). 



that the coupling strength is 51 = 0.3 x 2n MHz n(J; 
the difference in magnetic moments is A M sw 0.19/io |10| . 
where po is the Bohr magneton. The results of Fig. 
are for 1/B = 400/is/G 10], which corresponds to 
£ = (lnA0 2 A MJ B/(M7 2 ) w 7.9 for N = 10 5 atoms per 
species. Of course, even the reduced-space model [l3| 
cannot handle N = 10 5 atoms, but for £ = 7.9 then 
N = 10 2 will actually underestimate the N = 10 5 results. 
Hence the already good agreement between Fig-EJd) and 
Fig. ^ would actually improve if resources were available 
to manage the correct number of particles. 

We pause briefly to justify the ideal gas model. The 
collisional interaction strength is roughly A = 2Trhpa/m, 
where a is the off-resonant atomic s-wave scattering 
length. The 40 K scattering length is a = 176ao [13] > with 
ao the Bohr radius. For a typical density p ~ 10 13 cm~ 3 , 
the collisional coupling strength, in units of the atom- 
molecule coupling, is |A|/f2 sa 10~ 3 . Collisions should 
therefore be broadly negligible. In particular, a system 
of Fermi atoms coupled to Bose molecules is formally 
identical to a system of only bosons 0, and collisions 
are negligible for bosons under such conditions 0. 

Also, it should be noted that, because we have cho- 
sen Jl as the frequency scale, the above results are ac- 



tually for a resonance (atom-molecule coupling) of arbi- 
trary strength. However, the model © is broadly equiv- 
alent to the two-mode model in coherent association of 
condensate [13J, and it is well known that strong cou- 
pling can lead to dissociation to modes lying outside the 
two-mode system 0, 0], so-called rogue dissociation. 
Nevertheless, if the sweep is directed from above to be- 
low threshold, then rogue dissociation is negligible and 
the two-mode model is a good approximation. Hence, 
the above results are expectedly reasonable to describe a 
sweep across an arbitrarily strong resonance. 

Before closing, we turn to a related matter of principle: 
the nature of above-threshold molecules. Below thresh- 
old (d < 0), Fourier analysis delivers the binding en- 
ergy, huj B < 0, of the Bose-condensed molecules [lfj, 
lub — S — H'(ljb) + it] = 0, where £'(wb) is the finite 
self-energy of the Bose molecules and rj = + . Tun- 
ing the system above the two-body threshold (6 > 0) 
gives an imaginary u>b, and the bound state ceases to 
exist; nevertheless, Fig.|3shows a large AT = 1 molecular 
fraction. This result is not really a surprise, since the 
fraction of molecules must vary continuously from zero 
to unity across threshold. We conclude that any theory 
in which molecules abruptly cease to exist at threshold, 
while useful in their own right (e.g., for modeling bind- 
ing energies [Tol ] ) , are not a good rule of thumb for pre- 
dicting the existence of above-threshold molecules. Our 
interpretation is that, as usual in cooperative behavior, 
a macroscopic number of particles respond as a unit to 
a given external drive, thereby mimicking one- or two- 
body physics. At the least, this implies that many-body 
effects are sufficient but not necessary for the existence 
of above-threshold molecules. Moreover, we see in Fig. El 
that the above-threshold molecular fraction for d/Cl < 2 
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FIG. 3: Fraction of above-threshold molecular condensate for 
N — 100 (from Fig. |3Jc), dashed) compared to ground-state 
N = 1 (solid line) . Evidently, above-threshold molecules exist 
in the absence of many-body effects and, oddly enough, such 
effects can act to suppress the molecular fraction for positive 
detunings S/Q > 2. 
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is actually suppressed for the many-body case N = 10 2 . 

Of course, the idea of many-body stabilization of 
above-threshold molecules generally arises in the context 
of equilibrium thermodynamics, whereas the collisionlcss 
model describes non-equilibrium processes. The many- 
body suppression of the above-threshold molecular frac- 
tion may or may not carry over to the collisional regime 
(although we find elsewhere that, to a certain degree, it 
may |21|). However, we expect that the two-body equi- 
libration time is sufficiently long compared to the atom- 
molecule conversion timescale that, even in the presence 
of collisions, a two-body system can always be considered 
out of equilibrium; hence, the two-body ground state of 
Fig. should apply to the collisional regime as well. 

Conclusions. -We have investigated saturation in col- 
lisionless rapid adiabatic passage from a two-component 
degenerate Fermi gas to a Bose-Einstein condensate of 
molecules. Saturation indeed arises, but can be system- 
atically eliminated by introducing the timescale appropri- 
ate to cooperative interference effects. Physically, coop- 
erative interference effects arise from adding up the vari- 



ous pathways coupling the states having N — m molecules 
and m dissociated atom pairs with the states having one 
more (less) molecules and one less (more) dissociated 
pair, and the timescale for iV-particle interference turns 
out to be ~ InN/Cl 13] . It then makes perfect sense that 
cooperative (near-unit-efficient and macroscopic) rapid 
adiabatic passage will only occur over a timescale that 
is commensurate with constructive interference. Next 
we saw that our zero-temperature model agrees semi- 
quantitatively with our mean-field model llOf. indicating 
that temperature 01 an d pair correlations |14| are-as of 
yet-the main obstacles to collisionless cooperative con- 
version to molecules with near-unit efficiency. Finally, 
whereas studies of Feshbach resonances for both fcrmions 
and bosons have implicated many-body effects in the ex- 
istence of molecules above the two-body threshold for dis- 
sociation, we find that it is not necessary to invoke many- 
body effects to explain the existence of above-threshold 
molecules. 

Acknowledgements— One of us (O.D.) kindly thanks 
the Magnus Ehrnrooth Foundation for support. 



[1] W. C. Stwalley, Phys. Rev. Lett. 37, 1628 (1976); E. 

Tiesinga, B. J. Verhaar, and H. T. C. Stoof, Phys. Rev. 

A 47, 4114 (1993). [10 
[2] C. A. Regal, C. Ticknor, J. L. Bonn, and D. S. Jin, Na- [11 

ture (London) 424, 47 (2003). 
[3] K. E. Strecker, G. B. Partridge, and R. G. Hulet, Phys. [12 

Rev. Lett. 91, 080406 (2003); J. Cubizolles et ai, Phys. 

Rev. Lett. 91, 240401 (2003). 
[4] M. Greiner, C. A. Regal, and D. S. Jin, Nature (Lon- [13 

don) 426, 537 (2003); S. Jochim et ai, Science 302, [14 

2101 (2003); M. W. Zwierlein et al, Phys. Rev. Lett. 

91, 250401 (2003). [15 
[5] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 

92, 040403 (2004). 

[6] M. W. Zwierlein et al, Phys. Rev. Lett. 92, 120403 [16 
(2004). 

[7] C. Chin et ai, Science 305, 1128 (2004); J. Kinast et al, [17 
Phys. Rev. Lett. 92, 150402 (2004); T. Bourdel et al, 
Phys. Rev. Lett. 93, 050401 (2004). [18 

[8] M. Tinkham, Intro, to Superconductivity, (McGraw-Hill, [19 
New York, 1975). 

[9] E. Timmermans, K. Furuya, P. W. Milonni and A. K. [20 
Kerman, Phys. Lett. A 285, 228 (2001); M. Holland, S. [21 
J. J. M. F. Kokkelmans, M. L. Chiofalo, and R. Walser, 



Phys. Rev. Lett. 87, 120406 (2001); Y. Ohashi and A. 
Griffin, Phys. Rev. Lett. 89, 130402 (2002). 
J. Javanainen et al, Phys. Rev. Lett. 92, 200402 (2004). 
J. Chwedehczuk, K. Goral, T. Kohler, and P. S. Julienne, 
e-print cond-ma t/0409192 

0. Dannenberg, M. Mackie and K.-A. Suominen, Phys. 
Rev. Lett. 91, 210404 (2003); M. Mackie et al, Phys. 
Rev. A 69 053614 (2004). 

1. Tikhonenkov and A. Vardi, e-print cond-mat/0407424 
E. Pazy, A. Vardi, Y. B. Band, Phys. Rev. Lett. 93, 
120409 (2004). 

This means we will be able to fully account for the atoms 
being fermions without, per say, stacking them in trap 
levels, and that densities are understood to be uniform. 
J. Javanainen and M. Mackie, Phys. Rev. A 59, R3186 
(1999). 

A. Vardi, V. A. Yurovsky, and J. R. Anglin, Phys. Rev. 
A 64, 063611 (2001). 

A. Ishkhanyan et al, Phys. Rev. A 69, 043612 (2004). 

M. Mackie, K.-A. Suominen, and J. Javanainen, Phys. 

Rev. Lett. 89, 180403 (2002). 

J. L. Bohn, Phys. Rev. A 61, 053409 (2000). 

M. Mackie and J. Piilo (unpublished). 



